Analytical expressions for the local-field factor G(q) and the exchange-correlation 
kernel K xc (r) of the homogeneous electron gas. 
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We present an analytical expression for the local field factor G(q) of the homogeneous electron 
gas which reproduces recently published Quantum Monte-Carlo data by S. Moroni, D.M. Ceperley, 
and G. Senatore [Phys. Rev. Lett. 75, 689 (1995)], reflects the theoretically known asymptotic 
behaviours for both small and large q limits, and allows to express the exchange-correlation kernel 
K xc analytically in both the direct and reciprocal spaces. The last property is particularly useful in 
numerical applications to real solids. 
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The static local field factor G(q) of the homogeneous 
electron gas (HEG) is an important quantity because it 
represents the extent to which the particle interactions 
affect the exchange and correlation properties of this ide- 
alized system. The importance of G{q) in calculating the 
properties of real materials stems from the fact that it 
can be used as a key input in density functional calcu- 
lations fll§]. In fact, approximations of the unknown 
exchange-correlation energy functional of real (inhomo- 
geneous) systems involve the second functional derivative 
of E xc [n] , or exchange correlation kernel 

K xc (n ; |r-r'|) = / | K „, (1) 

where n is the HEG density. On the other hand, the 
local field factor and the exchange correlation kernel are 
simply related in Fourier space by ||: 

K™(q) ee J d 3 re-^K xc (r) (2) 

= -v c (q)G(q), 

where v c (q) — Aire 2 /q 2 is the Coulomb potential. 

Various approaches to obtain expressions of the local 
field factor were investigated in the past fi|-|l0||. Moroni 
et al. |4| obtained, by quantum Monte Carlo (QMC) sim- 
ulation, the zero temperature local field factor G(q) of the 
homogeneous electron gas at densities n Q corresponding 
to r s = 2, 5, 10 {n — 3 / (4na 3 B r 3 ) , where as is the Bohr 
radius). They also gave an analytical formula for G(q) 
which fits the QMC computed value, and has the right 
asymptotic limits at small and large wave vector, i.e.: 

G(q) ~ AO 2 for q - 0, (3) 

where Q = q/kp, hp is the Fermi wave vector, 
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fj, c being the correlation contribution to the chemical po- 
tential; and 
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G{q) ~CQ 2 +B 



for q -> 00, (5) 



where 

and e c is the correlation energy per particle. In our cal- 
culations we use the parametrization of B from Ref. [Q : 

(3 + b\x + b 2 x i ) 

1 /2 

where a; = r s , Oi = 2.15, a2 = 0.435, &i = 1.57, and 
b 2 = 0.409, valid for r s in the range 2-10 Q. 

The Q 2 behaviour of the local-field factor at large Q 
has been overlooked for a long time. It has been demon- 
strated in Ref. jl2| , and a clear discussion of its origin is 
presented in Ref. Q): its coefficient C is related to the 
change in kinetic energy in going from non-interacting 
(Kohn-Sham) electrons to interacting (real) electrons. 

In this paper we fit the values of G(q) computed in 
Ref. p[] in such a way to obtain a simple analytical ex- 
pression of both K^(q) and K xc (r), where K xc {r) is the 
exchange correlation kernel in real space. Our formula for 
G{q) is based on Lorentzian and Gaussian functions, and 
reads: 

G(q) = CQi + -^+aQ*e-™ 2 1 (8) 

g + Q 

where g = B/(A — C) and the two parameters a and 
P are fitted in order to minimize the differences with 
the numerical results of Ref. Q|. In particular, the best 
results are obtained taking: 

1-5 ^ 
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(3=^. (10) 
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Note that the Lorentzian contribution in Eq.^| is a sim- 
ple Hubbard-like term. This term alone yields already 
qualitative agreement with the numerical data of Ref. 
j|. Adding the Gaussian term allows us to reproduce 
quantitatively the numerical evaluation of Ref. g. In 
panels a, b, and c of Fig.l, G(q) given by Eq|8| is com- 
pared with the QMC results of Ref. Q for the unpolar- 
ized electron gas. The agreement is satisfactory, and is 
globally of the same quality as that obtained with the 
interpolation formula originally proposed in Ref. In 
panel d of Fig.l, we extend the comparison to QMC data 
for the fully spin polarized electron gas at r s = 100 p3| . 
Despite the fact that this value is well beyond the range 
considered in Ref. [Q for the parametrization of B (Eq.7), 
and that in Eq.8 we neglect spin polarization effects, the 
agreement is still fairly good. 
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Using Eq.|2|, and after Fourier transforming, we obtain 
the expression of the exchange correlation kernel K xc (r) 
in real space: 
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In Fig. 2, we compare this form of i"C xc (without the first 
term, which contains a three-dimensional delta function) 
with the (numerical) FT of the kernel derived from the 
Utsumi and Ichimaru parametrization of G(q) |^|. Be- 
sides the very desirable property of allowing passage from 
real to reciprocal space and vice-versa without numeri- 
cal transforms, the present form has another advantage, 
which is particularly useful in calculations for real solids: 
the absence of long-range oscillations, at variance with 
the case of those G(q) which contain a logarithmic sin- 
gularity for q = 2kp, as the Ichimaru- Utsumi one (see 
Fig. 2). This singularity is a peculiar property of the 
homogenous electron gas, originating from its spherical 
Fermi surface and from the absence of a gap between 
filled and empty states. Even in the HEG, its existence 
is not certain, and it is likely not present in real mate- 
rials flifl . Hence it is better, given the present level of 
ignorance about the exchange-correlations kernel of real 
materials and the computational difficulties arising from 
such a singularity, to use expressions of G{q) which do 
not contain it, and consequently do not yield the long 
range oscillations of K xc {r). 

In conclusion, we have presented a new parametriza- 
tion of published QMC results for the local field fac- 
tor G(q) of the homogeneous electron gas. Our analyt- 
ical form fits the numerical data with the same accu- 
racy as the form originally proposed in Ref. Q, has the 
right limiting behaviours for large and small q, and has 
the additional advantage of being analytically Fourier- 
transformable, a property which greatly simplifies its use 
in density- functional calculations for real materials. 

We are grateful to Saverio Moroni and Gaetano Sen- 
atore for providing us their unpublished QMC data for 
r s = 100, and for a critical reading of the manuscript. 
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FIG. 1. Panels a, b, and c: local field factor G(Q), 
(Q = g/fei?) as computed according to Eq. in compari- 
son with the QMC data for the unpolarized HEG of Ref. 
Panel d: Eq.8 compared with QMC data for the fully spin 
polarized HEG @. 

FIG. 2. The exchange-correlation kernel K xc of the homo- 
geneous electron gas as parametrized in the present work (full 
line), compared with the K xc of Ref. |B| (dotted line), in 
atomic units, plotted in reciprocal and direct space (upper 
and lower panel, respectively, with Q = q/kF and R = rfcj?) 



